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Abstract: We report about very high energy (VHE) 7-ray observations of the Crab Nebula with the 
MAGIC telescope. The 7-ray flux from the nebula was measured between 60 GeV and 9 TeV. The energy 
spectrum can be described with a curved power law gf = /o (S/300 GeV) (a+i,loglo(B/300 GoV)) with 
a flux normalization / of (6.0 ± 0.2 stat ) X 10~ 10 cm -2 s _1 TeV~\ a = -2.31 ± 0.06 sta t and b = 
-0.26 ± 0.07 st at- The position of the IC-peak is determined at (77 ± 47) GeV. Within the observation 
time and the experimental resolution of the telescope, the 7-ray emission is steady and pointlike. The 
emission's center of gravity coincides with the position of the pulsar. Pulsed 7-ray emission from the 
pulsar could not be detected. We constrain the cutoff energy of the spectrum to be less than ~ 30 GeV, 
assuming that the differential energy spectrum has an exponential cutoff. For a super-exponential shape, 
the cutoff energy can be as high as ~ 60 GeV. 



Introduction 

The Crab Nebula is the remnant of a supernova ex- 
plosion that occurred in 1054 A.D. at a distance of 
~ 2kpc. It is one of the best studied non-thermal 
celestial objects in almost all wavelength bands of 
the electromagnetic spectrum from 10~ 5 eV (ra- 
dio) to nearly 10 14 eV (7-rays). The radiation 
from radio to 7-rays (E < 1 GeV) is interpreted as 
synchrotron emission of relativistic electrons and 
positrons. At higher energies it is believed that 
inverse Compton scattering is the dominant gen- 
eration process of 7-rays |[8]|6l[T2). There is lit- 
tle doubt that the engine of the nebula is the pul- 
sar PSR B0531+21 (hereafter Crab pulsar), which 
is also a strong source of 7-rays detected up to 
10 GeV. 

In very high energy (VHE) 7-ray astronomy the 
Crab nebula was first detected with large signifi- 
cance at TeV energies by the pioneering Whipple 
telescope fl2l . Since then the Crab nebula was ex- 
tensively studied by ground based experiments at 



energies above a few hundred GeV. However, be- 
tween 10 GeV and ~ 200 GeV, observations are 
sparse. 

Here we present highlights of an 16 hour observa- 
tion of the Crab nebula and pulsar that was per- 
formed with the MAGIC telescope between Octo- 
ber 2005 and December 2005. A more detailed dis- 
cussion of the analysis and results presented here 
can be found in (3j. 

After a short description of the MAGIC telescope 
and the performed observations we present results 
from the analysis of the VHE-7-ray emission from 
the nebula and the search for pulsed emission from 
the pulsar. The paper is closed with some conclud- 
ing remarks. 

The MAGIC telescope 

The MAGIC (Major Atmospheric Gamma Imag- 
ing Cherenkov) telescope is located on the Canary 
Island La Palma (2200 m asl, 28.45°N, 17.54°W). 
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MAGIC is currently the largest single dish (17 m 
diameter) imaging atmospheric Cherenkov tele- 
scope (I ACT). The faint Cherenkov light flashes 
produced in air showers are recorded by a cam- 
era comprising 577 photomultiplier tubes (PMTs). 
The central PMT is modified for optical pulsar 
studies ifTOl . 

The current configuration of the MAGIC camera 
has a trigger region of 2.0 degrees in diameter 0. 
Presently, the trigger energy range spans from 50- 
60 GeV (at small zenith angles) up to tens of TeV. 



Analysis results of the Crab Nebula 
Source Morphology 

The morphology of the 7-ray emission was studied 
by generating sky-maps in three different bins of 
energy. The center of gravity (CoG) of the 7-ray 
emission was derived from the sky-maps by fitting 
them with a 2D-Gaussian of the form 



a ■ exp 



(x - x) 2 + (y - y) 2 



2a 2 



(1) 



where F ICS is introduced to account for a possible 
constant offset of the background subtracted sky- 
map. The CoGs obtained from the fitted x and y 
are shown in Figure Q] superimposed on the com- 
posite image of optical, IR and X-ray observations 
of the Crab nebula. The three measured CoGs are 
compatible among each other and coincide with 
the position of the pulsar. Note that the systematic 
uncertainty of the position is about 1'. 

The extension of the 7-ray emission region is com- 
patible with a point-like source. Upper limits on 
the source extension were calculated with a con- 
fidence level of 95%. The results for energies 
- 250 GeV (< 2.4') and > 500 GeV (< 1.6') are 
presented in Figure Q] In both cases the emission 
is constrained to originate from within the optical 
nebula. 

Energy spectrum 

Figure [2] shows the differential flux measurements 
by MAGIC multiplied by the energy squared, 
i.e. the spectral energy density distribution. A pa- 
rameterization of the spectrum with a power-law 



Figure 1 : The crosses indicate the CoG of the VHE 
7-ray emission at different energies (simple cross 
> 500 GeV; dot - 250 GeV; square - 160 GeV), 
overlaid on an optical image from HST and an X- 
Ray image from Chandra. The position of the Crab 
pulsar is marked with a black star. The error bars 
indicate the statistical uncertainty in the position 
of the CoG. The upper limits on the 39% con- 
tainment radius of the 7-ray emission region are 
indicated with circles (dashed > 500 GeV; solid 
- 250 GeV). 

ansatz results in a \ 2 of 24 for 8 degrees of free- 
dom. A better parameterization is obtained with a 
curved power-law ansatz. 

^ = f ( £ /3 ooGeV) (a+ " os - (£/300GcV)) 
dE 

(2) 

yielding a flux normalization /o of (6.0±0.2 sta t) x 
10~ 10 cm- 2 s- 1 TeV- 1 , a = -2.31 ± 0.06 sta t and 
b = -0.26 ± 0.07 s tat ± 0.2 syst . The X 2 of the fit 
is 8 for 7 degrees of freedom. 

For energies above 400 GeV the derived spectrum 
is in good agreement with measurements of other 
air Cherenkov telescopes. At energies < 400 GeV, 
below the threshold of previous measurements by 
IACTs, we compare our results with integral flux 
measurements obtained by [7] and ifTTIl . 

At lower energies one expects a continuous soften- 
ing of the spectrum with increasing energy. How- 
ever, this could not be demonstrated by earlier 
measurements. We derived spectral indices at ~ 
150 GeV, - 300 GeV, - ITeV and - 2.5 TeV 
from the flux measurements as well as from the 
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Figure 2: Spectral energy distribution of the 7-ray 
emission of Crab nebula. Below 10 GeV, measure- 
ments are by EGRET. In VHE 7-rays measure- 
ments are from ground-based experiments. The 
dashed line is a model prediction by (TJ. 

aforementioned results of the curved power-law fit. 
The results, shown in Figure [3] together with sev- 
eral predictions, indicate a clear softening of the 
spectrum with increasing energy. 

The predicted GeV 7-ray emission has a peak in 
the SED-representation (see Figure |2). If one as- 
sumes that the energy spectrum around the peak 
can be described with a curved power-law, we can 
determine the position of the peak from the curved 
power-law fit to be at 77 ± 47 sta t t« T 8yBt GeV. 
Within statistical uncertainties, the flux of 7-rays 
was constant over the entire observation period. 
Tested timescales were between a few minutes 
up to months. The average integral flux above 
200 GeV is (1.96 ± 0.05 sta t) X 1CT 10 cirr 2 s" 1 . 

VHE-7-ray emission from the pulsar 

Motivated by a possible pulsed 7-ray component 
at TeV energies |9), we searched for pulsed emis- 
sion in five bins of reconstructed energy between 
60 GeV and 9 TeV. However, a signature of peri- 
odicity was not found in any of the tested energy 
intervals. Derived 95% confidence level flux limits 
are shown in Figure |U 

We also performed a periodicity analysis opti- 
mized for a search of pulsed emission close to 



Figure 3: Spectral index derived from differen- 
tial flux measurements (dots) and from the curved 
power-law fit (black solid line, the dashed line 
gives the la confidence interval); Predictions by 
[1| (blue curve), [4| (green curve) and Q (red 
curve) are also shown. 

the threshold of the experiment (analysis thresh- 
old 60 GeV). Figure [5] shows the resulting pulse 
phase profile together with EGRET data above 
5 GeV and measurements in optical by MAGIC. 
An excess in VHE-7-rays is evident at the position 
of the inter-pulse in the same phase range where 
EGRET detected pulsed emission above 100 Me V 
(shaded region) and above 5 GeV. The MAGIC and 
EGRET > 5 GeV pulse phase profile match with a 
probability of 87%. One calculates a significance 
of 2.9 (X of the excess if the phase regions where 
EGRET detected pulsed emission above 100 Me V 
(shaded regions) are defined as signal region and 
the remaining phase intervals as background re- 
gion. 

The observed excess is not sufficient to claim the 
detection of a pulsed signal, therefore, upper limits 
on the number of excess events where calculated 
with a confidence level of 95%. Using the limit on 
the number of pulsed excess events we constrained 
the cutoff energy of the pulsar spectrum to be less 
than 27 GeV, under the assumption that the break 
in the energy spectrum can be described with an 
exponential cutoff. In case the energy spectrum is 
attenuated super-exponentially, cutoff energies up 
to ~ 60 GeV are allowed from our observations. 
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Figure 4: Upper limits on the pulsed gamma ray 
flux from the Crab pulsar; upper limits in differ- 
ential bins of energy are given by the blue points. 
The upper limit on the cutoff energy of the pulsed 
emission is indicated by the dashed line. 

Concluding remarks 

Here we reported on the currently most detailed 
study of the VHE 7-ray emission of the Crab neb- 
ula below 500 GeV. Most of the aforementioned 
studies were done in this energy region for the first 
time. Results from this study among others are a: 

• measurement of the differential energy spec- 
trum down to 60 GeV, which clearly deviates 
from a pure power-law behavior 



• determination of the inverse Compton peak 

1 1 

syst 



at77±47 stat +_ 1 4 ° 6 7 GeV 



• point-like emission region with a CoG coin- 
ciding with the position of the pulsar 

• constraint of the cutoff energy of the pulsar 
spectrum of < 27 GeV, assuming an expo- 
nential cutoff 
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Figure 5: Pulse phase profiles of the Crab pulsar. 
Lower figure: optical observations by MAGIC; 
middle figure: observations by EGRET above 
5 GeV; upper figure: pulse phase profile ob- 
tained by MAGIC. The shaded regions indicate the 
EGRET measured positions of the pulsed emission 
for 7-ray energies above 100 MeV. 
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